! ***************************************
FINISH
/CLEAR
/FILN, E53 
/TITLE, BUCKLING ANALYSIS OF ELLIPSOIDAL HEAD WITH VERTICAL SUPPORT LOAD
/UNITS,SI
! ********* 趨*********
DI=1000                    !Ͳ弰Բͷھ
T=10                       !Ͳ弰Բͷں
T1=18                      !ͷ
DI=1000+T                  !о
HC=DI                      !Ͳ峤
DR=630                     !֧װֱ
DS1=150                    !ֱ
D2=108                     !ֱֹ֧
NH=30                      !ʷ
NN=20                      !ͷʷ
FYY=4.27E6                 !ֵṹ趨ôֱغ
EXX1=2E5                   !ϵĵģ
NUXY1=0.3                  !ϵĲɱ
SY=235                     !
!****************ǰ***************************
/PREP7 
ET,1,181                   !趨Ԫ
R,1,T,T,T,T                !ǵԪʵ
R,2,T1,T1,T1,T1            !ǵԪʵ
MP,EX,1,EXX1               !ϵĵģ
MP,NUXY,1,NUXY1	           !ϵĲɱ

TB,BISO,1,1,2,  	   !˫Ըͬǿ
TBTEMP,0                   !¶ֵΪ0
TBDATA,,SY,0,,,,           !
!****************ģ***************************
CSWPLA,11,1,0.5,1,         !Բϵ
K,3,-DI/2,,,               !Բ˵
K,4,DI/2,,, 
L, 3,4                     !Բ
CSYS,0                     !ֱϵ
K,5,DI/2,-HC               !Ͳ˵
L,4,5                      !Ͳĸ
K,1,0,0,0                  !תĹؼ
K,2,0,HC,0                 !תĹؼ
L,1,2                      !
LPTN,1,3                   !߶λ
LDELE,5, 6, ,1             !ɾ߶
LDELE,4                    !ɾ߶εĹؼ
ALLS                       !ȫѡ
AROTAT,2,7 , , , , ,1, 6,60, ,	!תԲͷͲĸɷͷ
WPOFF,DR/2                 !ƶƽ
WPROT,0,90                 !תƽ
CYL4, , ,DS1/2, , , ,-DI   !ɸԲ

ASLV                       !ѡ
ASEL,R,LOC,Z,-DS1,0        !ѡ
ASEL,U,LOC,Y,0             !ȥ¶
ASEL,U,LOC,Y,DI            !ȥ϶
*GET,ACON,AREA,,NUM,MIN    !ȡи
ASEL,S,LOC,Y,0,DI          !ѡͷ
*GET,ACON1,AREA,,NUM,MIN   !ȡͷ

APTN,ACON1,ACON            !滥֣26
VDELE,1, , ,1              !ɾԲ

LSEL,S,LOC,Y,DI            !ѡи϶߶
ASLL                       !ѡ
ADELE,ALL, , ,1            !ɾѡ(9)
LSEL,S,LOC,Y,0             !ѡY=0λ߶
LSEL,R,RADIUS,,DS1/2       !ѡи¶߶
ASLL                       !ѡ
ADELE,ALL, , ,1            !ɾѡ棨10

CYL4, , ,D2/2, , , ,-DI    !ֹ֧ܸԲ
LSEL,S,RADIUS,,D2/2        !ѡ뾶ΪD2/2߶
ASLL                       !ѡ
ASEL,U,LOC,Z,0,D2/2        !ȥλô0D2/2 
*GET,ACON,AREA,,NUM,MIN    !ȡиţΪ5
ALLS
APTN,7,ACON                !ָ
VDELE,1,,,1                !ɾԲ
ADELE,10,11,,1             !ɾи

!**************************************************
LSEL,S,LOC,Y,-HC,0         !ѡλY=-HC0֮߶
LSEL,U,LOC,Y,-HC           !ȥλY=-HC߶
LSEL,U,LOC,Y,0             !ȥλY=0߶ΣѡLINEΪ12
LESIZE,ALL, , ,NH,, , , ,1 !趨Ͳ徭ʷ
LSEL,S,LOC,Y,-HC           !ѡλY=-HC߶
LSEL,A,LOC,Y,0             !ѡλY=0߶ΣѡLINEΪ45
LESIZE,ALL, , ,NN,, , , ,1 !軷ʷ
ALLS
LESIZE,3, , ,2*NN,, , , ,1 !߶ʷ
LESIZE,19, , ,NN,, , , ,1 
LESIZE,20, , ,NN,, , , ,1 
LESIZE,16, , ,NN,, , , ,1 
LESIZE,22, , ,3,, , , ,1
LESIZE,23, , ,3,, , , ,1  
LESIZE,18, , ,2*NN,, , , ,1
LCCAT,3,5                  !ε߶Ϊһ߶

WPROT,0,0,90               !תƽ棬Աʷ
ASBW,6                     !ʷֹ֧ڵ
AMESH,1                    !Ͳӳʷ
AMESH,8                    !ͷǵӳʷ
ASEL,S,,,9                 !ѡֹ֧
AATT,,2                    !
LCCAT,8,11                 !߶
MSHKEY,1                   !趨ӳʷ
AMESH,9                    !ӳʷ

ASEL,S,,,2,3,1             !ѡֹ֧ڵ
AATT,,2                    !
AMESH,2                    !ֹ֧ڵ
AMESH,3

ALLS                       !ȫѡ
ARSYM,Z,ALL, , , ,0,0      !Գ
NUMMRG,ALL, , , ,LOW       !ϲ
NUMCMP,ALL                 !ѹ
FINI

!*******************************************
/SOLU
ANTYPE,STATIC              !
OUTRES,,1                  !趨ѡ
NLGEOM,ON                  !η
PRED,ON                    !ԶԤ
KBC,0                      !Լ
NSUBST,20                  !趨Ӳ
ARCLEN,ON,5                !

LSEL,S,LOC,Y,-HC           !ѡͲ϶
DL,ALL, ,ALL,              !ȫԼ
LOCAL,11,1,,,,,90          !ֲϵ
LSEL,S,LOC,Y,60            !ѡṹԳλϵ
LSEL,A,LOC,Y,-60
LSEL,U,LCCA                !ȥ
DL,ALL,,SYMM               !ʩӶԳԼ
ALLS

LSEL,S,RADIUS,,D2/2,(D2+DS1)/4	!ѡ֧
NSLL,,1                         !ѡڵ
*GET,NTT,NODE,NUM,COUNT         !ȡѡڵ
NFYY=FYY/NTT                    !ֱغƽ
F,ALL,FY,-NFYY                  !ѡڵʩӴֱغ
ALLS
SOLVE
FINI
!*******************************************
/POST1

SET,LIST,2                      !бȷٽغ



